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It is shown that in ultrahigh energy inelastic neutrino-nucleon(nucleus) scattering the cross sec¬ 
tions for the boson-hadron(nucleus) reactions should exhibit geometric scaling on the single variable 
TA = Q 2 /Qs,t,A- The dependence on energy and atomic number of the charged/neutral current 
cross sections are encoded in the saturation momentum Q aa t,A. This fact allows an analytical com¬ 
putation of the neutrino scattering on nucleon/nucleus at high energies, providing a theoretical 
parameterization based on the scaling property. 
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I. INTRODUCTION 

One important property of the nonlinear perturbative 
QCD approaches for high energy deep inelastic ep(A) 
scattering is the prediction of the geometric scaling. 
Namely, the total 7 *p(A) cross section at large ener¬ 
gies is not a function of the two independent variables 
x and Q 2 , but is rather a function of the single variable 
t a = Q 2 /QLt,A- As usual, Q 2 is the photon virtual¬ 
ity and x the Bjorken variable. It was demonstrated Q 
that geometric scaling is the exact asymptotic solution 
of a general class of nonlinear evolution equations 111 
and it appears as a universal property of these kind of 
equations. The specific scaling solutions correspond to 
traveling wave solutions of those equations. The satu¬ 
ration momentum Q 2 at A (x; A) oc ~ A a x~ x 

(a ~ 1/3, A ~ 0.3) is connected with the phenomenon of 
gluon saturation. It is expected that the rise of the gluon 
distribution function at small values of Bjorken x be nat¬ 
urally tamed by saturation and circumvent the subse¬ 
quent violation of unitarity. In principle, geometric scal¬ 
ing is predicted to be present only on process dominated 
by low momenta. However, it was theoretically demon¬ 
strated 1 that the geometric scaling is preserved by the 
QCD evolution up to relatively large virtualities, within 
the kinematical window Ql at (x) ^ Q 2 ^ Qs at (x)/Aq CD . 
That is, the scaling property extends towards very large 
virtualities provided one stays in low- 2 ;. This kinemat¬ 
ical window is further enlarged due to the nuclear en¬ 
hancement of the saturation scale. In Ref. JJ, the ob¬ 
servation that the DESY-HERA ep collider data on the 
proton structure function F 2 are consistent with scal¬ 
ing at x < 0.01 and Q 2 < 400 GeV 2 was for the first 
time presented. Similar behavior was further observed 
on electron-nuclei processes Geometric scaling has 
been also observed in inclusive charm production [ 7 ] and 
predicted to be present also in heavy quark production 
in lepton-nuclei scattering 

Recently, the high energy lepton-hadron, proton- 
nucleus and nucleus-nucleus collisions have been related 
through geometric scaling ;9j. Within the color dipole 
picture and making use of a rescaling of the impact pa¬ 
rameter of the 7 *h cross section in terms of hadronic tar¬ 


get radius Rh , the nuclear dependence of the 7 *A cross 
section is absorbed in the A-dependence of the saturation 
scale via geometric scaling. The relation reads as 
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where the nuclear saturation scale was assumed to rise 
with the quotient of the transverse parton densities to the 
power A = 1/S. The nucleon saturation momentum is set 
to be Q 2 at (x) = {xq/x) x GeV 2 , where xq = 3.04 • 10" 4 , 
A = 0.288 and x = x \l + (4 m 2 /Q 2 )], with to/ = 0.14 
GeV, as taken from [l(j. The nuclear radius is given by 
Ra = (1.12A 1 / 3 — O. 86 A -1 / 3 ) fm. The following func¬ 
tional shape (scaling curve) for the photoabsortion cross 
section has been considered based on theoretical motiva¬ 
tions jf}: 


(riot ( t p ) = ^0 [7 e + r (0, (3) + In (/?)], 


P = a/Tp , (3) 


where 7 e is the Euler constant and T (0, fj) the incom¬ 
plete Gamma function. The parameters for the proton 
(nucleon) case were obtained from a fit to the small-x 
ep DESY-HERA data, producing a = 1.868, b = 0.746 
and the overall normalization was fixed by do = 40.56 
/tb. The parameters for the nuclear saturation scale 
were determined by fitting the available lepton-hadron 
data using the relation in Eq. 0 and the same scal¬ 
ing function, Eq. 0. One obtains <5 = 0.79 ± 0.02 and 
ttR 2 = 1.55 ± 0.02 fm 2 . The value of <5 favors a nu¬ 
clear dependence Q 2 at A oc A 4 / 9 Q 2 at , which it implies 
a faster increasing in comparison with the usual ansatz 
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The theoretical and phenomenological achievements 
summarized above have direct consequences on the com¬ 
putation of the cross sections of interaction of neutri¬ 
nos with nucleons and atomic nuclei. Ultrahigh energy 
(UHE) neutrino reactions are sensitive upon the behavior 
of the nucleon/nuclear structure functions at extremely 
low -x values and relatively large scales of electroweak bo¬ 
son virtualities. Namely, they can probe the domain 
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x ~ m\ vz /E l , ~ 10 8 at A;, ~ 10 12 GeV and virtu¬ 
alities Q 2 ~ ~ 10 4 GeV 2 , where mw,z are the 

boson masses. Therefore, accurate predictions for UHE 
neutrinos require precise extrapolation of the structure 
functions to the small-a:. Recently, the approaches con¬ 
taining saturation effects and/or scaling property have 
been compared and their main features were identified 
and investigated mill A very important feature is 
that, within the color dipole framework El , the charged 
(CC) and neutral (NC) current structure functions are 
described by the same mathematical expressions as the 
proton structure function up to a different coupling of 
the electroweak bosons. Therefore, the geometric scaling 
property should be present also in neutrino scattering 
on hadron targets and allows to obtain the dependences 
on energy and atomic number of CC/NC cross sections, 
which are encoded in the nuclear saturation momentum. 
In fact, geometric scaling holds for a large region of phase 
space and thus it contributes for an important part of the 
integrated cross section. Furthermore, the upper limit for 
scaling window should be enhanced by a factor ~ 10 — 100 
in the nuclear case. In what follows, the weak boson- 
hadron/nucleus cross section is computed and it is shown 
to exhibit geometric scaling on the scaling variable ta- 
The CC and NC structure function are further studied 
and compared with the available high energy approaches. 
The scaling property will allow to obtain simple analyt¬ 
ical expressions for the CC/NC neutrino cross sections 
at high energies. They will be used to provide theoret¬ 
ical parameterizations for the UHE neutrino scattering 
on nucleons and nuclei. 


II. UHE NEUTRINO CROSS SECTION 


Deep inelastic neutrino scattering can proceed via 
(charged current interactions - CC) or Z° (neu¬ 
tral current interactions - NC) exchanges. The stan¬ 
dard kinematical variables describing them are given by 
s = 2 mjsiEv (center-of-mass energy squared), Q 2 (boson 
virtuality), Bjorken x and y = Q 2 jxs (inelasticity vari¬ 
able). Here, itin is the nucleon mass and E v labels the 
neutrino energy. At small-a:, a successful framework de¬ 
scribing QCD interactions is provided by the color dipole 
formalism , which allows an all-twist computation of 
the structure functions. The physical picture of the in¬ 
teraction is the deep inelastic scattering (DIS) at low 
x viewed as the result of the interaction of a color qq 
dipoles, which are fluctuations of the electroweak gauge 
bosons, with the hadron target. The interaction is mod¬ 
eled via the dipole-target cross section, whereas the boson 
fluctuation in a color dipole is given by the corresponding 
wave function. The DIS structure functions for neutrino 
scattering read as EH, 

F? C l ' NG (x, Q 2 ) = ^ a tot (W ± (Z°) N - X ), (4) 



FIG. 1: Geometric scaling for the electroweak boson-nucleus 
cross section as a function of the scaling variable ta for dif¬ 
ferent atomic numbers (normalized to nucleon). 


crtot (W ± (Z°) N) = J d 2 r j dz \^ W± A° | 2 <7 d i P , (5) 


where r denotes the transverse size of the color dipole, z 
the longitudinal momentum fraction carried by a quark 
and V , 'r)f are the sum over dipoles of the wave func¬ 
tions of the charged or neutral gauge bosons, respec¬ 
tively. Their explicit expressions can be found for in¬ 
stance in Refs. El IT2I Il4l. Here one considers only 
four flavors ( u,d,s,c ) assumed to be massless. Heavy 
quarks (b, t) give relatively small contribution and will 
be disregarded. Color dipoles contributing to Cabibbo 
favored transitions are ud(du), cs (sc) for CC interac¬ 
tions and uu, dd, ss, cc for NC interactions. The dipole 
cross section adip (x, r; A), describing the dipole-target 
interaction, is substantially affected by saturation effects 
at dipole sizes r <1 1/Q sa t- 

Based on the fact that the expressions for the photon 
wavefunction and the electroweak gauge bosons, appear¬ 
ing on Eq. < 0 , are exactly the same up to the different 
coupling of the bosons to the quark color dipole, we can 
simply write, 
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where a em is the QED constant coupling, e f is the electric 
charge of the quark of flavor /. The constant A' qu i ra i = 
(A 2 + L 2 d + R 2 + R 2 d ) is the sum of the chiral couplings 
expressed as functions of the Weinberg angle 9w- 

The scaling present in the lepton-hadron cross section 
at high energies, as quantified by Eq. m and further ex¬ 
perimentally demonstrated, is automatically translated 
to the UHE neutrino scattering on nucleon or nucleus. In 
Fig. Qone presents the boson-hadron cross sections, Eqs. 
cap , as a function of the scaling variable ta for distinct 
nucleus as well as for the nucleon. They are normalized 
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FIG. 2: The charged (left panel) and neutral (right panel) 
current structure functions times R^/R 2 A versus Bjorken-a: for 
different atomic nuclei at Q 2 — m^y {m%). 


to the nucleon. The cross sections exhibit geometric scal¬ 
ing, verifying a transition in the behavior on t a of the 
cross section from a smooth dependence at small t a and 
an approximated 1/r behavior at large t a ■ Similarly to 
the lepton-hadron case, the transition point is placed at 
t a = 1. The asymptotic 1/r 6 dependence reflects the 
fact that the cross section scales as Qsat/Q 2 nrodulo log¬ 
arithmic corrections, with energy dependence driven by 
the saturation scale. The mild dependence at t a ^ 1 
corresponds to the fact that the cross section scales as 
oc oo log(Q 2 at /m 2 ) towards the photoproduction limit. 
The main features present in the cross sections, which 
are driven by the scaling function Eq. C9> can be qual¬ 
itatively reproduced in the phenomenological saturation 
models (see discussion in Ref. 0 )- The expressions get 
simplified to o-^ 011 oc r A 7 [1 + log (r A )] ( 7 ~ 1). 

In Fig. [3 one shows the charged and neutral current 
structure functions, given by Eq. 0 , rescaled by the 
factor Rp/R\ as a function of x and distinct nuclei, in¬ 
cluding the proton case. They are computed at scale 
Q 2 = (rn'z) and plotted down to extremely small 
x ~ 10~ 8 . For the proton case, the results are similar 
to the saturation model, whereas stay a factor about 1.7 
smaller than the unified BFKL-DGLAP approach with 
screening effects and a factor 3 in the case without screen¬ 
ing (for more detail on these results, see Refs. HI EH). 
This has consequences in the reduction of the magnitude 
of the neutrino cross section due to the nonlinear pertur¬ 
bative QCD effects. The saturation effects slow down the 
rise of the structure functions as x decreases. These ef¬ 
fects are sizeable even at larger values of order Q 2 ~ m^y 
because they also contribute to the leading twist part of 
the structure functions due to the geometric scaling win¬ 
dow. Concerning the nuclear structure functions, the nu¬ 
clear effects at the scale Q 2 ~ rri 2 v are relatively smaller 
than the usual expectation oc A 1//3 F^ V . As in the 
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FIG. 3: The neutrino-nucleus CC (left panel) and NC (right 
panel) cross sections as a function of neutrino energy E v and 
for different atomic nuclei. The curves are rescaled by the 
factor R 2 /R A . 


present calculation Q 2 at A oc A 4 / 9 , a smaller nuclear 
shadowing should be expected. Therefore, this implies 
a smaller reduction in the magnitude of the neutrino- 
nucleus inelastic scattering due to nuclear shadowing in 
contrast with recent approaches using the strength of the 
nonlinear term proportional to A 1 / 3 . 

The total CC (NC) neutrino-nucleon cross sections as 
a function of the neutrino energy and atomic number are 
given by the integration over available phase space and 
read as, 
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where a minimum <5^ in oc 0(1) GeV 2 is introduced in 
order to stay in the DIS region and Gp is the Fermi con¬ 
stant. Here, one considers UHE neutrinos, where the 
valence quark contribution stays constant and physics is 
driven by sea quark contributions. Hence, the xF^ C ’ NC 
contribution should be negligible and it will be disre¬ 
garded. In Eq. 0 , the nuclear dependence on the struc¬ 
ture functions is implicit. 

In Fig. [3 one shows CC and NC neutrino cross sec¬ 
tions, given by Eqs. 0, rescaled by the factor R 2 /R 2 A 
as a function of neutrino energy E v and distinct nuclei, 
including the nucleon. Based on the color dipole picture, 
one used the simple relation Fp « (2/ll)F 2 in Eq. 0, 
which gives a very small contribution. It is worth men¬ 
tioning that the calculations are also shown for either 
low neutrino energies ( E v ^ 10 11 GeV), where a part of 
the contribution to the integrated cross section would be 
out the scaling window. For the nucleon case, the results 
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where A/.,) are overall normalizations, C[ l \ are numerical 
constants with i = CC, NC, w sca i = 6 A and a = b/S. 
This implies in a mild power-like rise w SC ai — 0.2 for the 
neutrino cross section in contrast with other theoretical 
approaches. The nuclear dependence is approximately 
linear, cr™f lel oc Acr™ cleon , once b ~ 6 and hence a « 1. 
The remaining constants are given by, 


A 7 ' i) Hr 


o'oGpTn^ vz \ f a Xn 


8tt 3 X 


B 


W 


O.Q 


E 


p 2 1 

f e f 


(ii) 


Ci = 7r(2mjv)" scal (rri^z) 


CSC(7T!/ Sca l) (1 - l/ sca l) , 


= 7T (rriW' Z ) b CSC (nb) (1 - b ), 


FIG. 4: Comparison with other approaches: unified DGLAP- 
BFKL, NLO DGLAP and CGC calculations (see text). 
Present result (geometric scaling) given by the solid line . 


are close to the saturation models 0 . For the nuclear 
case, the results indicate a weak nuclear shadowing, as 
discussed above. A comparison with other high energy 
approaches is presented in Fig. 0 for the nucleon case 
(see Ref. m for detail on those results). At energies 
E v <; 10 12 GeV, the current calculation produces a re¬ 
duction by a factor 2 in relation to both NLO DGLAP 
and unified BFKL-DGLAP approaches and a milder en¬ 
ergy dependence. On the other hand, the result has sim¬ 
ilar behavior as the CGC model. The reason is that the 
CGC model takes the geometric scaling property and its 
extrapolation to the saturation region in the dipole cross 
section. 

The property of the geometric scaling allows to ana¬ 
lytically calculate the UHE neutrino cross section given 
by Eq. ©. The integration on Bjorken-r; can be car¬ 
ried out by making a change of variables. The integra¬ 
tion on x is then replaced by an integration on ta of 
the scaling function, Eq. This produces the explicit 
expression da/dQ 2 oc (3 3 F 3 [ 1 ,1, 1 ; 2, 2, 2; —/?]£”“, where 

Anin = Q / Qs&t = Q /^)i An ax = Q /Qsati^ = 1) and 
3 F 3 is a hypergeometric function. The further integra¬ 
tion on Q 2 can be performed, using the remaining leading 
term on virtuality and supposing y = Q 2 /xs -C 1. The 
result using the leading term represents with good accu¬ 
racy the full numerical calculation. The ratio between 
the full calculation and the leading term is independent 
of energy and takes a constant value for E„ ^ 10 8 GeV, 
giving R COI = (T;,(exact)/cr„(approx.) = 0.82. Therefore, 
the analytical calculation provides a theoretical param¬ 
eterization for the UHE neutrino-nucleus cross sections 
based on the geometric scaling property. They are as 
follows, 
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where one uses the notation Bcc = 4, Bnc = Achiral 
and zz scal = b — w sca i. Numerically, this gives a total cross 
section = 1.48 x 10 “ 34 A“(£; y /GeV ) 0 - 227 cm 2 . The 
cross section above can have implications for neutrino 
observatories because experiments are planned to detect 
UHE by observation of the nearly horizontal air showers 
in Earth coming from neutrino-air interactions [11] • A re¬ 
duced cross section produces a smaller event rate for such 
neutrino-induced showers and could compromise the de¬ 
tection signal. However, the rate of up-going air showers 
initiated by muon and tau leptons produced in neutrino- 
nucleon reactions just below the surface would increase, 
being possibly larger than the horizontal air shower rate. 

In summary, it is demonstrated that the cross sections 
for FU ± (^°)-nucleus processes in UHE inelastic neutrino 
scattering should exhibit geometric scaling on the sin¬ 
gle scaling variable ta = Q 2 /Q 2 at A - This implies that 
such dimensionless scale absorbs their energy and atomic 
number dependences. Based on this property, an analyt¬ 
ical calculation at high energies is made possible. This 
allows to propose a theoretical parameterization for the 
UHE neutrino cross sections supposed to be valid at ener¬ 
gies E v 10 8 GeV. The resulting cross section has a mild 
energy behavior on energy in comparison with the usual 
QCD calculations based on linear evolution equations. 
Moreover, the nuclear dependence is approximately lin¬ 
ear on the atomic number. 
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